Mesoscopic motion of atomic ions in magnetic fields 
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We introduce a semiclassical model for moving highly excited atomic ions in a magnetic field 
which allows us to describe the mixing of the Landau orbitals of the center of mass in terms of the 
electronic excitation and magnetic field. The extent of quantum energy flow in the ion is investigated 
and a crossover from localization to derealization with increasing center of mass energy is detected. 
It turns out that our model of the moving ion in a magnetic field is closely connected to models for 
transport in disordered finite-size wires. 

Interacting particle systems in strong magnetic fields show a rich variety of complex phenomena. The source 
of this complexity is competition between the magnetic and Coulomb interactions which are of inherently different 
character. With changing strength of the external field the corresponding systems undergo a metamorphosis involving 
qualitatively different states. In atomic physics attention focused for more than a decade on the hydrogen atom in a 
magnetic field (see ref. Q] and references therein) on which detailed experimental and theoretical investigations yielded 
many beautiful insights into semiclassical and quantum aspects of nonintegrable systems and significantly enhanced 
our understanding of the new features arising due to the presence of the external field. With increasing degree of 
excitation and/or increasing field strength the electronic motion of the classical atom shows a transition from regular 
to irregular, i.e. chaotic, behavior and intermittency. More recently it has become evident that the non-separability 
O i' of the collective, i.e., center of mass (CM), and electronic motion of atoms in the presence of a magnetic field leads 
to a variety of two-body phenomena. The corresponding coupling of the CM and electronic motion is fundamentally 
different for neutral and charged systems. Examples of two-body effects in neutral systems are the classical chaotic 
diffusion of the CM Q or the existence of weakly bound giant dipole states || . For atomic ions the interaction of the 
CM and electronic degrees of freedom is more intricate and manifests itself in a continuous, classical flow of energy 
£SJ ' from the collective to the internal motion and vice versa. Detailed studies of the classical dynamics of rapidly moving 
highly excited iJe + -ions in a magnetic field showed that this energy exchange leads to the self-ionization process ||| of 
the ion. Very little, however, is known about the quantum properties and behavior of moving highly excited atomic 
ions. 

Using Landau orbitals for the CM motion in zeroth-order and fixed nucleus zero-field wave functions for the 
electronic motion to estimate their coupling matrix elements, it was demonstrated || that there exist a number of 
different physical situations for which the interaction between the collective and electronic motion becomes strong. 
The latter induces a strong mixing of the CM and electronic motion and is a potential source of interesting quantum 
properties of the ion, in particular when the dynamics of the corresponding classical ion is chaotic. A detailed 
investigation of the coupled CM and electronic motion of the highly-excited quantum mechanical ion in this regime, 
O v which is the subject of interest of the present paper, is however a highly nontrivial task: we are dealing with five 
nonseparable and strongly mixing degrees of freedom in a regime of very high level density which depends on a number 
. \ of parameters (field strength, total energy, etc.). The ab initio description of the quantum dynamics in the above 
regime goes even beyond modern computational possibilities and we thus seek a model approach that captures the 
essential physics of the problem. Here we propose and analyze a semiclassical model of the excited ion, and explore 
the consequences of coupling between its CM and electronic degrees of freedom. 

Since we deal with the interaction of the CM and electronic motion in atomic ions we first have to introduce 
collective (CM) and relative variables in the Hamiltonian describing the atom. The total pseudomomentum K |(| is 
a conserved quantity associated with the CM motion which, in spite of the fact that its components perpendicular 
to the magnetic field are not independent, i.e. do not commute, can be used to transform the Hamiltonian to a 
particularly simple and physically appealing form Q which for the iJe + -ion reads TL = Tii + Ti.2 + H3 where 
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where to, Mo and M are the electron, nuclear and total mass, respectively, a = (Mq + 2m) /M and Q is the net 
charge of the ion. B is the magnetic field vector which is assumed to point along the z-axis. (R, P) and (r,p) are 
the canonical pairs for the CM and internal motion, respectively. H involves five degrees of freedom since parallel to 
the magnetic field the CM undergoes free translational motion and can be separated completely. 

Hi and H3 depend exclusively on the CM and electronic degrees of freedom, respectively. Hi describes the free 
motion of a CM pseudoparticle with charge Q and mass M. H3 describes the electronic motion in the presence of 
paramagnetic, diamagnetic as well as Coulomb interactions which, in analogy to the hydrogen atom [jlj, exhibits 
an enormous complexity of classical and quantum properties with changing parameters, i.e. energy and/or field 
strength. Hi contains the coupling between the CM and electronic motion of the ion and represents a Stark term 

with a rapidly oscillating electric field jj ^B x — ^B x Rj^j determined by the dynamics of the ion. It is the 

interaction Hamiltonian H2 which is responsible for the interesting quantum effects that will be investigated in the 
present Letter. 

^From the above it is natural to consider the representation of the coupling Hamiltonian H2 in a basis which consists 
of products of eigenstates $cm of Hi and \& of H3. Calculating the corresponding matrix elements we encounter 
some selection rules which are of immediate relevance to our model (see below). Since the total angular momentum 
component parallel to the magnetic field C z is a conserved quantity for H and since the corresponding CM angular 
momentum Lcai, and electronic angular momentum L z are conserved quantities for Hi and H3, respectively, the ma- 
trix elements of Hi involve only CM and electronic states with magnetic quantum numbers which are correspondingly 

different. In addition we have the relation < ®' CM \ \P - %B x Rj \§cm >= iM(E' - E) < $' CM \R\$ C m > which, 

together with the dipole selection rules for electronic transitions, allows only changes of the CM (/1) and electronic 
(to) magnetic quantum numbers by one and requires a change of energy for the CM motion. According to ref. |g] the 
matrix elements of Hi involve a factor y/N for N » \fi\, (N is the Landau principal quantum number of the CM 
motion) which yields the scaling of the coupling matrix elements of Hi with respect to varying CM energy. 

Our model for the moving H e + -ion in a magnetic field is built up from three key constituents associated with the 
Hamiltonians Hi, Hi and H3. The equidistant and infinitely degenerate spectrum of Hi is completely characterized 
by the CM quantum numbers N and \x. Hi represents the integrable part of the system which is coupled via Hi 

to the chaotic part represented by H3. The classical dynamics of H3 |?j depends on the scaled energy and angular 
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momentum which are given by L z = L z IjB) 3 and E = E (2B) 3 , respectively. In order to locate the regime of 
chaotic electronic motion we have made Poincare surfaces of section (PSOS) of the classical dynamics of H3 (^E,L z ^j 

for a dense grid of values of the scaled energy and angular momentum. Starting with a completely chaotic phase 
space for L z — we find that with increasing values of the angular momentum the fraction of chaotic phase space 
volume decreases rapidly. In contrast, the regime of negative L z that yields predominantly chaotic phase space is 
much larger. If we take a typical scaled energy of E = —0.1, for example, and require more than 90% of phase space to 
be occupied by chaotic trajectories we obtain the regime L z = [—2.71, 0.136] which corresponds to L z = [—200, +10] 
for B = 10~ 5 (we use atomic units throughout, i.e. the field strength B — 1 a.u. corresponds to 2.35 • 10 5 Tesla). The 
mixing of electronic eigenfunctions belonging to negative values of the angular momentum L z represents therefore an 
important 'open channel' with respect to the coupling of the chaotic electronic motion to the CM motion. In this 
region of L z the CM motion of the classical ion is strongly affected by coupling to the chaotic electron. Indeed, studies 
of the classical dynamics of the ion close to the ionization threshold demonstrated 0] that large negative values of the 
angular momentum are an inherent feature of intermittent dynamics as well as a prerequisite for the self-ionization 
process, whereby energy transfer from the CM to the electron results in ionization. We therefore turn our attention 
to this subspace and investigate the quantum mechanical properties of this channel. 

The spectrum belonging to the chaotic Hamiltonian H3 will be represented by a random matrix ensemble, which 
is the appropriate semiclassical description ]s|Jl0[|. Since the Hamiltonian H3 possesses a generalized time reversal 
invariance ||, which consists of a rotation by n around the i-axis and a subsequent conventional time reversal 
operation, the proper ensemble is the Gaussian Orthogonal Ensemble (GOE). While the GOE provides the fluctuations 
of the chaotic levels, we still need to specify the mean level density (MLD) as a function of the energy, field strength 
and in particular the angular momentum L z . Our approach to the MLD is via the semiclassical Thomas-Fermi 
formula. Performing the appropriate scale transformations, fixing L z , transforming to cylindrical coordinates and 
subsequently performing the integrations over <p,p z ,p p , z we arrive at the following result for the semiclassical MLD 
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with A = ( jt- ( + -^z ) — -E 1 ) — p 2 ■ The remaining integration over p has to be performed numerically. Starting 
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with L 2 = a general feature of the MLD is its rapid decrease with increasing L z whereas for negative values of 
L z its decay is much weaker. Not only the fraction of chaotic phase space but also the absolute phase space volume 
persists down to large negative values of the angular momentum, to L z w —250, and therefore pl_ (E, B) represents in 
this regime the density of irregular states. Having specified our GOE, whose MLD at the center of the band is given 
by eq.(2) providing the levels of H3, we turn to the calculation of the coupling matrix elements introduced by 7i 2 - 
The size of the matrix elements of H 2 can be determined from a semiclassical relation between off-diagonal matrix 
elements of an operator, and the Fourier transform of its classical autocorrelation function [ jTT| . The variance ct^ 
of the matrix elements of H 2 depends on the energies of the states they couple and is very small when the energy 
difference is greater than the level spacing of Hi For states close in energy, cr^ 2 appears as 
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where selection rules determine the numerical coefficient of eq.(3) [|l0|. Efficient evaluation of the above equation 
through classical trajectory calculations is by no means trivial and is done by introducing a suitable ensemble average 
via the periodogram approximation fl^| . Fig. 1 shows x = PlJch 2 \i as a function of the angular momentum, 
calculated using eqs.(l-3). As examples we plot ./V=250, 300 and 400 (at m = 0) and p < 0. Due to the selection 
rules N increases (decreases) by 1 as m decreases (increases) by 1. As seen in the figure, the largest values of x ue 
mainly in the interval m as [—20, —200]. 

Based on the structure of H the organization of our model is as follows: We have an array of sites each of which 
corresponds to eigenstates of Hi and H3 labeled by particular values of N, m and p. We thus assign to each site 
the CM quantum numbers N, p of a Landau orbital belonging to Hi, and the levels of one member of the GOE 
representing the irregular levels of H3, which are labeled by m. The energy levels at each site are the H3 levels plus 
the level of Hi- They are coupled by random matrix elements to levels of their neighboring sites, as imposed by the 
selection rules. The variance of the random elements <r^ 2 is given by eq.(3). The model thus has a one-dimensional 
structure, where sites labeled by N, m and p are comprised of levels taken to be from a GOE whose density is given 
by eq.(2). The ion model so defined resembles the semiclassical pump model of Arnold diffusion fl3(| . The classical 
stochastic pump model describes Arnold diffusion in terms of pumping of otherwise regular trajectories via weak, 
irregular motion within the Hamiltonian system JlJ]. Similarly, the classical CM motion of the ion is coupled to the 
chaotic motion of the electron due to the magnetic field. The semiclassical pump model of Arnold diffusion was found 
to be equivalent to models of single-particle transport in disordered wires, which predict localization |p^|-|l7||. The 
localization length of the semiclassical pump model was expressed in terms of corresponding classical parameters. 
Likewise, we should be able to predict any localization of the quantum ion in terms of the corresponding classical 
system. 

The values of all the ion parameters (N, p, E3, L z , B) are embodied in the model parameter x = Pl z Wh 2 1 1 examples 
of which are shown in Fig. 1. Though x clearly fluctuates over the sites of the model, we could nevertheless estimate 
the localization length for this model in terms of an average, x, over x- As seen in Fig. 1, the range in m over which 
energy transfer can most appreciably occur is m w [—20, —200], though the model also encompasses a wider range of 
m, where classically the ion is still predominately chaotic. In terms of the average x, the localization length, £00, of 
the model assuming an infinite number of sites, is £oo = 47r 2 % 2 [|l3Ul6| . Because the ion model has a finite length of 
about 200 sites, the localization length can be estimated using finite-size scaling arguments for band random matrices, 
like the ion model, as £ w £oo/(l + c£,oo/L), c w 1 where L is the length, or number of sites of the model. The 
N dependence of £ arises from ^ ~ N, since \<th 2 \ ~ V^/V (see above). Solving for pl z and \<jh 2 \ using eqs.(2) and 
(3) we find, e.g., that £ w £oo w 30 for A^=400 (m—0). Thus starting in a Landau level near N sa 400, quantum 
flow is restricted to about 30 sites in m. This is in contrast to the classical ion, where for a corresponding initial CM 
energy of m 6 x 10~ 7 a.u. there is no such restriction in the chaotic motion over L z . Since x varies as yN the ion 
remains localized, in contrast to the classical ion, for TV up to values near 4000. Taking the cyclotron frequency to be 
f .4 x f0~ 9 a.u., this corresponds to a CM energy of about 6 x f0~ 6 a.u. 

These arguments, while using well known predictions for equivalent random matrix models, nevertheless depend on 
our being able to use an average of x over the length of the model to estimate £. As a check, we have studied numerically 
the ion model to compare predictions of £ using the actual x which varies as determined by the semiclassical results 
for the ion calculated with eqs.(l-3), with results using %, defined as the average x over the length of the model. Our 
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numerical model ranges from m=-230 to 0. Each site consists of 12 levels of the GOE. Though this is a very small 
number, it is all we could include computationally while incorporating also the largest possible number of sites. We 
average over 6 realizations of each set of parameters and calculate the localization length of the eigenvectors. Results 
are shown in Fig. 2, where we see that both models, the ion and the simpler version with give the same results for 
£, and are close to the line £ = 0.72^00, where ^ = 4-7r 2 x 2 , and the factor 0.72 is purely an artifact of using only 12 
levels per site, which is seen upon comparing with numerical results where more levels per site were used |13| . 

The occurence of the crossover from localized to delocalized CM motion that we find for the ion at a CM energy 
of about 6 x 10~ 6 a.u. with the above parameters will of course vary with the strength of the field and the internal 
energy of the ion, as well as its mass, since we have here considered only He + . This crossover should be observable 
spectroscopically since the regime of very weak and strong mixing show inherently different level spacings and absorb- 
tion features. Finally, we mention that the above-investigated quantum mixing of collective and electronic motion for 
atomic ions is certainly of interest also for charged molecular systems in a magnetic field. Here the heavy vibrational 
and rotational degrees of freedom couple, for heteronuclear systems, to the collective motion of the molecular ion 
providing a potential source of new rotational and vibrational structures. 

This work was completed during a visit of P.S. at the Max-Planck Institute for Physics of Complex Systems which 
is gratefully acknowledged for its hospitality. DML acknowledges support from NSF CHE 95-30680. 



Figure 1. The model parameter \ = PL z \cTi. 2 \i where pr, z and \<7h 2 \ are calculated using eqs.(2) and (3), respectively, 
as a function of the angular momentum quantum number, m. At m — 0, TV = 250 (short dashes), 300 (long dashes) 
and 400. 

Figure 2. Localization length computed for the ion at different values of the Landau level, N, at m=0, from which 
X is computed with eqs.(2) and (3). Filled circles are results for the ion model and open circles for simpler version 
described in text. The line is a fit through the data. 
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